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SUMI4ABT 



In this paper solutions are presented for the determination 
of transient thermal stress distributions, in hollow circular 
cylinder* of elastic material, in which the temperature is a 
function of the radial coordinate as well as time. Indicated 
solutions are given for an externally insulated cylinder where the 
surface temperature is a known function of time, and for the case 
where the temperature-time relations are known for a heating fluid 
within the cylinder. Numerical results have been included for the 
specific case of a hollow steel cylinder subjected to the flow of 
saturated steam. ?or tho particular set of conditions considered 
the tangential and axial stresses were found to be a maximum at 
the heated surface at approximately five seconds after the start 
of heating and attained a magnitude of 52,500 pounds per square 
inch. Stresses of this magnitude if often repeated could cause 
fatigue failure and result in rupture of a line carrying any 
appreciable pressure. It is conceivable that for more extreme 
conditions the surface stresses might reach magnitudes high enough 
to cause a permanent set in the material and residual stresses 
would be present when the cylinder was again cooled to the initial 
conditions. Stresses of these magnitudes warrant consideration 
in the design of steam lines that may be subjected to repeated 
heating cycles. 
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IKTHQWCTIOS 



The thermal stress distribution in hollow circular cylinder* 
caused by a ter.pers.tnre gradient in the radial direction only 1* 
well known, The case of steady state heat flow in such cylinder* 
is treated in several textbook*. The problem of etrcss distri- 
bution for the case in which the temperature is a function of both 
the radial coordinate and tine lend* it*olf to solution by the sane 
stress equations, providing that the teraperature-tima relations are 
known. Analytic solutions for transient temperature distributions 
are given by many authors for the case where, fron an initial uni- 
fora tenperature, one or both of the surfaces are subjected to 
linear tenperature gradients, or are suddenly raised to a given 
value and thus maintained. The process for solution of s\ich prob- 
lens is to substitute into the already available stress equations, 
expressions for the tenperature a* a function of tine and radiu* 
for the given condition of heat flow. It can be seen then that 
for the case of the hollow circular cylinder, as for other confi- 
gurations, the general thermal stress problem is quite straightfor- 
ward and evolves itself, for the more complex conditions of heat 
flow, into that of the determination of tenperature distribution* 
as a function of time and distance* 

Consider now the more general problem of the thermal stress 
distribution in a cylinder subjected to the flow of a medium 
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through the cylinder different in temperature to that of the 
cylinder* In this case the inner and outer surface temperatures 
will vary with time and in general will not be linear. The solu- 
tion of such a problem also involves the consideration of a surface 
heat transfer coefficient which itself may have & magnitude which 
varies with temperature and consequently with time. Analytic solu- 
tions for the transient tenperaturo distribution is a cumbersome 
double series of Bessel functions which does not readily lend itself 
to numerical results. However, analytic solutions are available 
for unidirectional heat flow in the semi-infinite solid in terms 
of exponentials and the probability function, which are readily 
adopted to numerical results and ean be aoplied to the hollow 
circular cylinder within certain limitations. Such solutions ars 
valid for infinitely long thin cylinders and give a close approxi- 
mation to the temperature distribution for thick cylinders, in whieh 
ths wall thickness is not large compared to the radius. The use- 
fulness and practicability of such a method to determine transient 
thermal stresses in hollow cylinders when subjected to a heating 
medium is ths subject of this paper. 

In the development of the problem, two separate methods of 
attack were considered and their results compared. First, a numer- 
ical method for the transient temperature distribution based on a 
difference equation rather than the applicable differential equation 
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vas used to determine the time-temperature gradient* Secondly* 
a purely analytic method wae used to determine the time- tempera- 
ture relations and in both instances a numerical integrating process 
was used to evaluate the integrals in the stress equations. It 
might be noted at this point that the stress equations will give 
exact results, within the accuracy of the numerical integration* 
if the time-temperature relations thus determined are the true 
temperature gradients for the heat flow considered. 

It is the intent of this paper to indicato a practical method 
of solution for such thermal stress problems rather than to obtain 
accurate numerical results for any specific problem, however, in 
order to more clearly indicate the feasibility of this approach* 
a simple numerical problem has been carried through, using average 
values of thermal and olastic properties with slide-rule accuracy. 

The solutions of more complex problems are indicated as the method 
is developed. The problem chosen for its simplicity is that of a 
hypothetical ten-inch I. D. by twelve-inch 0. D. steam line, ini- 
tially at 100 °F suddenly subjected to a saturated steam flow having 
a temperature of 500°y* It is assumed that the Bteam starts to 
condense immediately and remains at the satixration temperature 
throughout the transient heat transfer process. The surface heat 
transfer coefficient was taken &e an average value over the inner 
surface of the pipe, considering film type condensation, and assuming 
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it not effected "by temperature, -he magnitude of the surface 
heat transfer coefficient was chosen from average experimental 
values for film type condensation as f oun& in the literature. The 
reason for thee -a "basic assumptions and their validity are further 
discussed in the "body of the text. 
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GEHEBAL I Z AT 1 015 OF PROBLEM AND BASIC ASSUMPTIONS 



The general differential equation of transient heat flow la 
obtained from the familiar "basic Fourier equation for the conduc- 
tion of heat ( <^ - — A^d.'t) and a heat "balance. Disregarding 
heat sources or sinks and equating the heat gain to the difference 
between the heat entering and heat leaving an elemental volume, 
the general differential equation expressed in rectangular coordin- 
ates is: 



b t c P e 



+ £ ("H?) 



If k, the thermal conductivity, is considered inverable with 
temperature change and the material is homogeneous and isotropic, 
k will be the same in all directions and constant. 




In general, engineering materials can be considered isotropic 
and homogeneous end this is an ordinary assumption made in engi- 
neering calculations. However, the assumption that K is constant 
may not be valid for all materials within the range of accuracy 
desired in the resulting calculations. Tho order of magnitude of 
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the variation of k with temperature my be seen by citing mild 
steal as an eraanle, which has a thermal conductivity of approxi- 
mately 26 BTtj/hr ft at 32 °T and 21 BTU/hr ft *F at 1100°F. 

The desired tiae-tenparafcure-dlatanco relations for the heat- 
ing or cooling of various shapes are obtained by the integration 
of the general conduction equation, substituting in the necessary 
boundary conditions for the case in question. One of the simplest 
cases is the unidirectional heat conduction in a plate of semi- 
infinite thickness, with the free surface at x * 0. Tf the 
term k/cpP is defined aa the thermal diffusirity and indicated by 
the symbol X the general conduction equation for unidirectional 
heat flow reduce* to} 



To solve this equation, it it necessary to determine the 
formulae which will ezprerts the initial ar.d boundary condition# 
which the temperature satisfies. If a semi-infinite eolid with 
its plane face at x s o t initially at the constant temperature of 
zero and suddenly subjected to a beating medium, which has allied 
with it a surface heat transfer coefficient E, the initial and 
boundary conditions are: 




Initial condition T= O Wrtfc-N t = o 
Boundary condition — 3JC +■ _tL (Ts - Tcu - 



o At x =0 
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or defining H/k &a h the boundary condition isx 

- + -K (Ts ~ To.') ~ o At X-O 

dX 

The solution of the differential equation for unidirectional 
boat conduction with these initial and boundary conditions will give 
the temperature as a function of time and distune* for the serai- 
infinite solid. Such a solution will also he valid for a plat* 
of finite thickness having a large ratio of surface to thickness, 
if it i* heated, as before, at the face x s * and the other face 
is perfectly insulated. To use such a solution for a hollow cir- 
cular cylinder, externally insulated, and subjected to the flow 
of a heating medium it nust be aseuraed that the curved section 
of the cylinder can under certain circumstances be considered as 
a flat insulated plate. This assumption will necessarily introduce 
an error in the transient temperature distribution thus determined 
for the cylinder and will be a function of the ratio of thickness 
to radius. This fact must be considered in the calculation of 
stress distributions for stresses thus deternined are only as reli- 
able as the teaperature distributions need in their calculation. 

To verify the validity of this assumption the applicable 
difference equations may bo set up for the plane plate and the 
hollow circular cylinder and the results compared. As explained 
under the heading "Numerical Solution of Temperature Distribution", 
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the msserical method converges rapidly to the analytic solution# 
fills bolng the ease* a comparison of numerical methods should also 
justify the uso of the analytic solution for the plane seal-infinite 
solid for stress calculations for the hollow circular cylinders 

within the Basse range of accuracy# 

Taking a section of the plane solid of unit width* AX thick* 
and dividing it up aa shown in the figure below* a hoat balance 
gives: 



Jk(T,-Tz)At + M-nj-Tx)At = (Tz-Ti) 

WHErRfc: *T^‘ = Tz at (t+At) 

Oe (r.-TO + (T-s-Tz) =: rl (T 2 '-T 2 ) 

WHeee. M = CpPax 7 - 

ThetivJ: T* - Ta + IT -zTa +T 3 
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A similar heat "balance for the circular section where the 



area "between section* progreasively increase* gives s 




Awa© 



FIGUI33 B 



Jk(T>-T x )&t +Jfe(i.o + A*Ad)(T*-T a )At = C P PA^(i +A]ca©)(T Zy l-T l) 

Vs/Me-RE-: T-J = I2 AF-T&R (t+at') lM CyUMDfcR 

Or (x-H) + (to 4AXA©HT5-Ti) - m (» + Ajxa ©) (Hr -lz) 



It can he seen that equation* 1 and Z are identical in fora, 
the only difference "being the addition of the term Axa© for the 
circular cylinder. For the ordinary problem, in order to satisfy 
the basic assumptions of such a numerical solution, ax must be 
chosen as a low order term. This means that the radius of the 
cylinder will be large as compared to ax and thus Ae must also 
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toe a low order term. The magnitude of these values may readily 
toe seen toy taking as an example a hollow cylinder of five inch 
inside radius and a ax 3 ,01 ft (Hote that this is the largest 
value of ax chosen for numerical calculation) thens 

A© = Tana© - oJL*i2 - .024 Bad. 

5 

and ax a© - .00024 (Which is essentially a 

second order tons). 

Then for the use of a modulus M ■ 4s 

Tj.r s + ' T> - ( 2 .ooo 24 )T\ + (». 00024^^ 

4. 0004 0 

A comparison of equations 3 and 1 shows that this point* 
namely at a radius of five inches* T^ and are practically 
identical. However, in the numerical solution for temperature 
distribution, the temperature at any point is used to determine 
the temperature at the next adjacent point, etc., through the 
solid. This means that the error in using equation 1 instead of 
equation 2 will toe an accumulative error and places a restriction 
upon the thickness of the cylinder for a given error. This fact 

perhaps can toe seen more clearly toy noting from Figure A th&tj 

/ 

= A* + hvAXA© 

Wheres n is the number of steps to toe taken 
for the given thickness chosen. 
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It should also he noted that the error is not a function 



of Ax alone hut a function of the product of ax and A© „ 
which essentially means that for a given radius cylinder and wall 
thickness the error involved will he the same regardless ef the 
choice of ax. With these facts In mind, it is apparent that the 
limiting factor is the ratio of wall thickness to radius regardless 
of their absolute magnitudes. 

Another pertinent observation is that at the first instant 
of heating the two equations almost coincide and the error in 
temperature distributions increases with each step. As the tempera- 
ture distributions thus determined, only enter the stress equations 



the error involved in stress calculations would be to plot tempera- 
ture times radius vs. radius curves for both method# and compare 
their areas. 

It is the author’s opinion that the error involved due to this 
basic assumption will be of the same order of magnitude or less 
than the errors inherently involved in thermal and elastic proper- 
ties, where ratios of wall thickness to radius found in ordinary 
steam lines are to be used. 




the true criterion for 
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NUMERICAL SOLUTION fOH T3»fEERAfU12E DISTRIBUTION 



Unidirectional heat flow will occur in a flat slab at a 
distance from tho edges and also in an insulated rod of uniform 
section* As far as temperature distribution is concerned, the size 
and shape of the rod is immaterial and is only of consequence in 
calculating total heat flow. Consider such a rod of unit width, 
a thickness ax and divided up into sections ax long as indicated 
in the figure below: 




If now a finite surface heat transfer coefficient H is con- 
sidered between the surface (0) and the ambient fluid temperature 
(a) then for the point (0) a heat balance gives: 

H (TV- To) At = - -fe (T ( -To) at + C P ?/yx.‘ (To - To) 
Vhere To* - To At (t -tcX) 
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The finite difference e<joation thus obtained is based 



upon the following assumptions: 

a. That At be taken aaall enough *o that there is a neg- 
ligible error in heat flow calculated by the tern 

(Tj-TgJ/nx. 

b. That A t be taken email enough so that there is a neg- 
ligible error in neglecting any but adjacent areas 
(ax by ax), 

c. That ax be taken small enough so that the teepsrature taken 
at the center of the element can be used to cocroute the 
heat capacity. 

Define s M = CpP A x 2 = ax 2 

At. K At 

M - HaY 



Using these relations: 



Nl(Ta-To) + (T -To) r (To -To) 



Solving for T,5 

To = 2jvJ.Ta + (\ - 2 H -4 2 \ To 4- 

k\ r\ r\ 

This equation nay be written in terns of "influence factors* 

as follows: , 

To - ^uoTU + p 00 To 4- f-o.t; 
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The subscripts indicate in order 0 the influence that the 
present temperature has on the point whose temperature at (t-VAt) 
is to he found for the period of time At. In other words, if the 
new temperature at point (0) is to he found after the time At, 

I M Indicates the Influence of point (a) on the temperature at 
point (0) during the time interval At. 

If is to depend upon ? 0 in a positive sense, then; 

Poo = l- 2M+2 ^ O 

Fd 

OR rd-2N + 2 2 O 

id ^ 2NI + 2. 

Since all of the quantities entering into K are inherently 
positive, it is apparent that M must he larger than two , for the 
case where a surface transfer coefficient is involved. Depending 
upon the magnitude of the surface transfer coefficient and the ther- 
mal constants of the material, this may place a limitation upon & 
convenient averaging process. The value of M^4- was found suitable 
for numerical temperature calculations for steel. 

The general method is to use the weighting process, as described 
above, for surface temperature calculations and the simpler avera- 
ging process for interior temperature determinations. It is con- 
venient to choose a scale such that fractions can be neglected. 
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If the fraction 1/2 was consistently added or dropped# an error 
would tend to accumulate. Therefore, a rule of dropping the l/2 
after even digits and adding 1/2 after odd digits can he adopted# 
as; 1/2-0# 3/2-2, 5/2-2, 7/2-4, etc. 

The choice of the modulus M is discussed in reference (5). 

The larger the modulus M the sooner the numerical results will 
converge to the analytic solution (i.e. in leas steps) and the 
greater the number of calculations necessary to cover the same time 
period. After due consideration of the work involved# a modulus 
was chosen to obtain rapid convergence to the analytic solution 
during the early time intervale. 

Now consider the heating of a one- inch thick plats Insulated 
on one surface and subjected to a heating medium on the other sur- 
face. For simplified calculation# consider I and B as constant. 

By choosing Beveral values of ax (and the compatible at) for a 
modulus M=4, surface temperature vs. time curves can. be plotted to 
determine the validity of the choica of the magnitude of ax. As 
the accuracy of the solution is based on the assumption that ax 
be email, such a plot should indicate any variation in the results# 
due to an increase in the magnitude of ax. 

For ordinary Steels 
I = ,50 ft 2 /hr 
k = 25 BTO/hr, ft, °F 
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Assume 



3 1000 'BTTJ/hr, ft 2 , °F (A fair approxim- 

tion for f il«~typ* condensing steaa — sse Appendix I) 

Sample determination of H» ^ao* ^oo> ^ol a ^ * 

Let AX = .01 F-T. 

Then: Nl = Hax 

M. 

=. AX x 
K At 

At = (.01^ _ 0.5 x icr* 

4K.5 

At - O-S a lO’ 4 "* 3boo - -18 5e-c. 

Fao - 1H — 2X0.4 - 0.2 

t^\ 4 . 

f*oo - (1 - 2N + 2 ) - (>~ O. S-4-2-0 ^ - 0.3 
r'v 4 

t“>o = - — — 3z- ~ O.S 
r'l 4 

Sidilarly for the other value* of Ax: 



AV -s 12 X .01 - .12 iNCHE-s 

— iooo_* .01 _ .4 o 
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Or At tr AX*. Wh&RC- K\=A 
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The equation* for the surface temperature "become! 



(1) 


To' 


= o.2 To. + 


o.3o To 


4- 


0.5 T, 


(2) 


To' 


- O.i 5 T<k + 


0.35 To 


4- 


0.5- T, 


0) 


To 


— O.loTo, 4- 


0.4.0 To 


4* 


o.«5T 


(4) 


To 


= 0.0 57a. + 


0.45 To 


4- 


0.5T, 



On set of calculations will he given at thiB point to indi- 
cate the form and procedure. The results only for other numerical 
calculations are included in Appendix II and plotted on the enclosed 
curve sheets. 
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AML mo SOIiTJTIOS FOR TZMPEJULTtJBE DISTRIBUTIOH 



To determine the tirae-teraper&ture-di stance relation* in the 
serai-infinite solid for the case where the initial temperature is 
aero and the surface temperature is a function of time <Mt)» the 
following equation is valid! 



OO 





© 



Also for the case of the serai-infinite solid, initially at 
sero temperature » which is heated at the surface x-3 by radiation 
froia a medlura at constant temperature T & , the applicable equation 
is: 




er-fc. * 



hx + h 1 fct 

_ e 



e.rfc 






Where erfc is defined as one sinus the probability function 
or error integral. 



e< fc * 



= ( 
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It can readily be shown by direct subs ti tut ion into the 
differential equation of unidirectional heat flow that these 
equations do satisfy the differential equation and initial and 
boundary conditions. 

Equation (l) was given to emphasise the fact that whenever 
the surface temperature of the semi-inf inite ,is known as a function 
of time cj>(t), and this function is continuous and can be integrated 
over the time interval (in equation (l)) to be considered* the 
resulting time-teaperature-distance relations for the solid are 
then determined. This opens ths avenue for the solution of a number 
of interesting problems. As notsd before, a numerical solution 
for surface temperature is readily adopted to thermal properties 
of the material which vary with temperature, as well ae a surface 
heat transfer that is a function of temperature. This method can 
also be extended to the oase where the temperature of the heating 
medium is some known function of tine, however, it should be noted 
that such a solution for surface temperature, under the above stated 
conditions may become very tedious, as the influence coefficients 
may have to be calctilatnd for each step of the process as well as 
choosing a At that is small. Although the numerical solution 
itself, if carried through the solid to the desired depth, will 
give the time-teisperatnre-distance relations directly, the desired 
surface temperature $(t) may be determined over a shorter interval 
with the resulting saving of much numerical work. 
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Figures (C) and (£>) show the determination of the surface 
teaper&ture relations 4>(t) for a short interval of time, Figure 
(£>) shows that the surface temperature vs* time curves very closely 
approximate a straight line on a log-log plot, and can readily he 
put into an analytic expression, 

Equation (2) has "been used to determine the time~tesiperafcur®« 
distance relations in a perf ectly- insulated, one-inch-thick plane 
solid of large surface area. The results of this can he used as 
a close approximation of the true transient temperature gradient 
In an insulated 10-ineh inside diameter hollow cylinder of one-inch 
wall thickness, 

for the ultimate calculation of the stress distribution, it 
is convenient to plot temperature— distance relations for various 
intervals of tine, in order to perform a graphic or numerical inte- 
gration of the integrals in the stress equation, Figure (6) shows 
a plot of these relations and shows how closely the numerical method 
parallels the analytic solution. Figure (£■) gives the equivalent 
tsaq'.erature tine* distance vs, distance relations, and indicates & 
graphic integrating process. It night he preferable and possibly 
a little more accurate to use a numerical integrating scheme such 
as the trapezoidal or Simpson's rule. 

Numerical results ere not as tedious as a first inspection of 
equation (2) might indicate. Tabulated values of the error integral 
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are available and many of the parameters are constant for one 



instant of tine. A simple tabulation of the quantities given 
helov will materially facilitate calculation, all other para- 
meters being constant for a given instant oi time. 

(1) V in Fr 



( 2 ) 

(3) 

W 

(5) 

( 6 ) 

( 7 ) 

( 8 ) 



TC 

\r*t 

e/fc TrT 

(-y|=r + ^t) 

h* 

h% + h 1 Kt 

e 

T 



The results of calculations made for the seal-infinite solid 
for a penetration of one inch are given in appendix to correlate 
the data used to plot curves. 
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BSTSEMINAT I C$T OP STE3SS DISTEIBUTIGN 



If the transient temperature distributions determined by either 
numerical or analytic methods for the seal-infinite solid under the 
conditions specified are considered good approximation* of true 
temperature distributions in the hollow cylinder, calculation of 
stress distributions are straight forward. At any instant of time 
the stress becomes a function of the thermal and elastic properties 
of ths material for the geometric conf iguration considered, and can 



be readily calculated. The following thermal stress equations for 
ths finite hollow cylinder are given by Timoshenko: 



These equations are based on thermal and elastic properties 
that do not vary with temperature. Within this limitation they 
are exact for an infinitely long cylinder of any wall thidmess 
and any given temperature distribution so long as it is in the radial 
direotion only. The principal planes of stress are In the radial. 
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axial, and tangential directions with normal stress acting on 
each and stress acting along a plane perpendicular to this axis 
are independent of the position of this plane. 

In order that the solution be exact when a cylinder of finite 
length be chosen , boundary conditions require that external forces 
equal in value and distribution to the stresses in the axial direction 
act on the ends. In the usual case such forces will not be present* 
so that the stresses near the ends will be very different than those 
calculated by this method, this effect diminishes as the distance 
from the ends is increased and practically disappears at a point 
sufficiently far away. Writers differ as to the distance from the 
ends where this paint is located and it actually depends upon the 
ratio of tho thickness to cylinder diameter. For the purpose of 
oar problem, this influence will be considered negligible. 

It is interesting to note that in the case where the surface 
of a solid, initially at uniform temperature, is suddenly raised to 
a given temperature and thus maintained, the maximum stress occurs 
at the surface at that instant. In fact, the maximum stress (i.a. 
tangential and axial stresses which are equal at the surface) will 
always occur at the heated surface for a solid initially at uniform 
temperature. The stresses thus produced axe equivalent to the stress 
produced in a thin plate for the ease temperature rise in which a 
change of area is prevented. This stress is only momentary and 
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rapidly drops a* the increase of temperature spreads inward and 
allows an expansion of the outer layer. Host published solutions 
have this common boundary condition, namely, an instantaneous 
surface heat source. 

For the case of radiation into a solid from a fluid, where 
the surface temperature is a function of time and equal to the 
initial solid temperature at the firBt instant of heating, the 
conditions are vastly different. The maximum surface stresses 
will not occur at the first instant of heating buit at some later 
time, and depending upon the conditions of heat flow may even occur 
at the steady state condition. For the case of the sudden rise of 
the surface temperature of a solid, initially at constant tempera- 
ture and ultimately cooing to a steady state condition, it is essen- 
tially going from an unstressed condition to a stressed condition. 

In the problem chosen, however, the cylinder at an initial uniform 
ternnorature with the outer surface perfectly insulated will, upon 
heating of the inner surface, be taken from one unstressed condi- 
tion to another as the cylinder finally reaches the temperature 
of the heating medium. It is reasonable to assume then that maxi- 
mum stress conditions will occur at some intermediate time. 

Calculated values of inner surface stress are tabulated in 
Appendix II for a one-inch thick, ten-inch I. D. circular steel 
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cylinder along with the stress distribution across it® section at 
the time of aaxiaum stress, The cylinder being insulated on its 
outer surface a initially at 100°P and suddenly subjected to a satur- 
ated steam flow which is maintained at 500°F. The results of these 
calculations are graphically represented by figures ( B ) » (^) and (O) 
and will be further discussed in the results and conclusions. 
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APPESDIX 2 



Accurate thermal stress determinations in steam line* are 
usually linked with a knowledge of the surface heat transfer co- 
efficient. The matter of the physical phenomena involved as well 
as the magnitude of sucl/ coefficients is the subject of a good 
deal of conjecture. This coefficient is defined by the heat flow 
equations 

q. = HA(T g -T a ) 

\sheres q. - Heat Plow 

H - Surface transfer ooefficient 
A - Transfer area 
T fc » Fluid temperature 
T»- Solid s.irface temperature 

Values of H determined by experimental means are in general 
more reliable than those determined analytically* and particularly 
for condensing steam. 

There are two types of condensation* namely, dropwl** and 
film condensation. Either or a combination of the two type* may 
exist, depending largely upon the cleanliness of the fluid, clean- 
liness of the condensing surface and the surface conditions of the 
material concerned In the condensation process. In any case, heat 
transfer coefficients for the dropwise process are In general many 
times higher than for film condensation, but are more difficult 
to determine. 
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Author# differ widely as to the magnitude of the difference between 
the two types* giving coefficient# from 4 to £0 or more time* larger 
for the none efficient dropwise process. It has been determined, 
that In general clean steam with or without non-aandansib.le ga«*» 
always condenses In a film on clean surfaces* whether rough or 
smooth. 

Experimentally determined values of heat transfer coefficients 
are given in terms of an average value over a prescribed area* 
discounting the fact that the value of the coefficient may vary 
considerably over the area depending upon the film thickness, Eusselt 
in an analytic investigation of heat transfer coefficients on cyl- 
indrical surfaces found the film thickness and consequently the 
transfer coefficient to vary around a horizontal tube. Hie analysis 
shows the minimum film thickness and consequent maximum transfer 
coefficient to exist a few degrees on either side of the top of a 
horizontal tube. Els investigation is equally valid for condensa- 
tion inside or outside of a tube. By integration around the tube 
an average coefficient was determined. Experimental determinations 
have shovn the average coefficient thus determined to be somewhat 
in error and correction factor curves are available for use with 
TTuaaelt *s equation to give representative values of the surface heat 
transfer coefficient for condensing steam on tubes. 
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If now dropvise condensation v^re to be considered, the 
stress problem immediately becomes very involved for tho local 
transfer coefficient would vary appreciably from droplet to un- 
wetted surface and tho heat flow would ne longer he unidirectional. 
For such a problem an entirely different approach would have to be 
considered. Tho value of H-1000 BTU/hr ft 2 used in the numerical 
calculations is an arbitrary value chosen from the references for 
film-type condensation of steam on a horizontal tube. 
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RESULTS iVD C'>'CLUSI^S 



It 1# apparent fro® th* magnitude of the Inner surface stress 
determined, for th® numerical problem considered, that thermal 
•tresses under serera Conditions of repeated heating cycles could 
contribute to ultimate rupture of a line under press-ire* The tan- 
gential stress on th® inner surface due to proa sure In the line 
will be tonsil® and to sots© degree counteract thermal stresses* 
However, for the 10~inch I. D. by I2«inah 0. D. line considered, 
the stress du® to a pressure of 515 f/in 2 in th© line would b® ap- 
proximately <f r m 500 #/ln 2 In compression and <fe * 1,640 #/in 2 in 
tension. These stresses are extremely snail in Magnitude as com- 
pnrod to the thermal stresses determined. 

It should be remembered that this paper was motivated by a desire 
to fornul&te a practical approach to the actual thermal stress 
problems in hollow circular cylinders which are extremely difficult 
if not impossible by parel. analytic approach. Kecesserily, many 
at sumptions were made, and to some degree they have been Justified 
in the tart of the paper. It is obvious then that the results of 
the numerical calculations should be substantiated by an experlaen- 
tal investigation before they are accepted as accurate enough for 
practical use* 

As little data is available on the Magnitude or variation ef 
surface transfer coefficients, fer condensing steam in hollow 
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circular cylinder* under high temperatures end pressure#, Uie 
value chosen could he considerably in error* as it is a value 
extrapolated fro® experimental results for condenser tubes under 
lover ter.pers.ture 8 and pressures. Xt is conceivable that a 
consideration o' the change of the magnitude of the surface trans- 
fer coefficient with temperature difference across the film as 
well as ths change in thermal properties with solid, temperature 
would materially affect ths thermal stresses determined. The 
simplifying assunptlcn that the line wae perfectly Insulated w^uld 
also constitute sons error in stress calculation* beyond th« time 
at which the outer surface temperature begins to rise. Although 
under steady stats conditions insulations aro available that are 
extremely effective, under transient conditions, temperature dis- 
tributions will necessarily be affected by th© heat capacity ef 
the insulation. This condition was not considered in the numeri- 
cal problem which injects soae error in tho resulting stress deter* 
isinationo for times beyond approximately five seconds. Perhaps, 
with some of tho qualifying assumption# in mind, this paper will 
stimulate further inv o s t igat ions along ths east ft lines toward a 
usable practical solution for thermal stresses in steam lines. 
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Further study of the practical, probien of thermal stress©* 

In hollow circular cylinders during actual transient boundary con- 
ditions falls Into the several broad categories listed below; 

(a) The experimental determination of ths true boundary 
conditions that «fiy eri s t in the actual problem during the 
transient state* such as those conditions that cay exist 
when a steam line is suddenly subjected to a flow of high 
temperature steam, 

(b) Jin experimental varif lent ion of the numerical example 
given in the text as to the tine at which maximum stress 
occurs* its magnitude, as well as the stress distribution 
across ths section at tints of maxima stress, 

(c) An experimental study of the magnitude of the surface 
transfer coefficient under high temperature* and pressures 
in hollow circular cylinders. How the value of such a co- 
efficient varies with pressure , tesmeratura difference 
across the condensate film, and angular position In a hori-< 
aontal cylinder. 

(d) A further justification of the us* of temperature dis- 
tributions as determined for unidirectional heat flow in 
the stress calculations for hollow circular cylinders. 
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Thors are two inunodiato avenues of approach to the last 
category, first, by direct comparison of the temperature dis tri- 
hut ion* as determined by merely nua rical methods fer the plane 
plate and hollow circular cylinder, and, secondly, such a tempera- 
turn comparison could be nade by the us© of a conformal napping 
process froa a plane finite aoction to a circular section* i. plot 
could then be nade of the error in stress calculations as & fvtac- 
tion of the ratio of wall thickness to radius, which would verify 
the usefulness of the approximats method for a problem of given 
geosnetry. 
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APPINDIX II 



tabulated values uaod to plot curvet 

Figure (A' - (Numerical method - Surface temp. ) 
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Figure (A'. - (Analytic method - Surface temp.) 
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Figure (E>^ - (numerical use t hod - Solid temp.) 
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Figure (e>) - Analytic Method - Solid Tenp 
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n«are (6) • (Analytic Solution - Surface street vs. tiae) 
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Figure (/=) - (Analytic Solution - Stress distribution at tine (5*0 sec) 
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